Generalized vorticity in transitional quantum turbulence 
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Transition to condensate state in the degenerate gas of bosons is studied using the Gross-Pitaevskii 
equation. It is shown that adiabatic invariance of an enstrophy (mean squared vorticity) based on a 
generalized vorticity c?irZ(|'i/)|^v)/|'i/)|^, controls dynamics of the transitional turbulence just before 
formation of the condensate with its tangles of quantized vortices (IV'Pv is the weighted velocity 
field defined on the macroscopic wave function -0). Scaling of the angle- averaged occupation number 
spectrum, A^^ ~ , has been obtained for the free developing transitional turbulence for the weak 
nonlinear (and completely disordered) initial conditions. Results of the three-dimensional numerical 
simulations have been used to support the theoretical consideration. 

PACS numbers: 03.75.Kk, 05.45.-a, 67.40.Vs 



Introduction. In this paper the term quantum tur- 
bulence is related to the stochastic dynamics described 
by the nonlinear Schrodinger equation for a macroscopic 
wave function that specifies the index of the coherent 
state (see, for instance, [H-[l3l)- Each coherent state 
evolves (to the leading order) along its classical trajec- 
tory. The Gross-Pitaevskii equation for instance, is 
used for description of such process (in a dimcnsionlcss 
form) 



^~dt 



(1) 



In particular, this equation gives a description of the for- 
mation of Bose-Einstein condensates from the strongly 
degenerate gas of weakly interacting bosons. 

It is believed that the nonlinear Schrodinger equation 
can also describe development of the coherent regime 
at a certain stage of evolution. Theoreti- 
cal considerations (see for instance, [l|-[ll) show that 
this regime can appear in a long-wavelength region of 
wavenumber space after the breakdown of the regime of 
weak turbulence: the Kagan-Svistunov (KS-) scenario. 
This is supported by recent tree-dimensional numerical 
simulations for the weakly interacting Bose gas (see for 
instance [l|). It is expected that there exist three differ- 
ent stages in the time development for the KS-scenario: 
weak turbulence, transitional turbulence and condensate. 
A pre-condensate (coherent regime) begins its formation 
in a long-wavelength region at the transitional turbulent 
regime. This transitional turbulent regime is the most 
difficult for theoretical description [3] , [1| , [HI . There is 
a problem to define a continuous vorticity filed in the 
quantum turbulence due to non-rotational nature of the 
velocity field defined on the wave function given by the 
Schrodinger equation. Velocity field defined on the wave 
functions is a potential field. Vorticity of such field van- 
ishes everywhere in a single-connected region. It is be- 
lieved that in the condensate state itself all rotational 
flow is carried by quantized vortices (the circulation of 
the velocity around the core of such vortex is quantized). 
This idea turned out to be very fruitful, and recent exper- 
iments give direct support for existence of such vortices 
in the condensate [IBl i [13] i [3 ■ However, for the transi- 



tional turbulence apparent absence of a continuous field 
associated with vorticity is a difficult problem. More- 
over, it is clear that the only flux of the particles from 
the region of lager energies toward the future condensate 
in not sufficient for transformation of the completely dis- 
ordered initial weak turbulence into the condensate with 
its local superfluid order and with the tangles of quan- 
tized vortices. A certain vorticity-like quantity have to 
be involved in this process just before appearance of the 
tangles of quantized vortices. 

In present paper we study a generalized vorticity 
defined on the weighted velocity field. The generalized 
vorticity is not vanishing in the bulk of the flow. Then, 
we show that adiabatic invariance of an enstrophy (mean 
squared generalized vorticity) controls the transitional 
turbulence just before formation of the condensate with 
its tangles of quantized vortices. 

First of all, let us recall certain properties of velocity 
fleld defined on the wave function. In order to calculate 
the expectation value of velocity one can evaluate the 
time derivative of space displacement (r) : 



(v) 



d(r) 
dt 



Substituting provided by the Eq. (1) one obtains 



(2) 



2m 



(3) 



where we have integrated by parts with corresponding 
vanishing or periodic boundary conditions. The inte- 
grand in the right-hand-side of Eq. (3) is the reaZ- valued 
distribution of the quantum velocity 



(4) 



The real-valued quantum velocity filed itself one can 
obtain comparing the definition (v) — J\r\^p\'^dx. with 
Eq. (3), that results in 



2m 



(5) 
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i.e. j = IfApv is the flux of probability or weighted veloc- 
ity. Representation (5) can be readily transformed into 



2m 



■0*. 



(6) 



i.e. the quantum velocity (6) is a potential field and 
curl V = everywhere in a single-connected region. 

Generalized vorticity. Let us define generalized vortic- 
ity on the weighted velocity field as 

curl (|?/'|^v) , Vl?/;!^ x v 
Lu = = curl V H '-^r^ (7) 



First term in the right-hand-side of Eq. (7) is the ordi- 
nary vorticity. When there is no the quantized vortices 
this term is zero due to Eq. (6). In this situation the sec- 
ond term in the right-had-side of Eq. (7) determines the 
generalized vorticity. In this case w is orthogonal to the 
"plane of motion" , which is stretched over the vectors: v 
and VjV'P, at any point of the space. In this sense the 
transitional turbulence can be considered as a locally two 
dimensional one. 

Dynamical equation for the generalized vorticity simi- 
larly to the dynamic equation for the density 



(8) 



does not contain the nonlinear terms explicitly 



dt 



= [Vip X V(AV'*)] + [VV-* X V(AV')] (9) 



(all the right-hand-side terms in the Eqs. (8), (9) come 
from the linear term of the Eq. (1)). This cancellation 
of the nonlinear terms in the dynamical equation for the 
density, (8), results (after integration by parts) in the 
conservation law for the total number of the particles: 
N — J \ip\'^dx = const. The cancellation of the nonlinear 
terms in the dynamical equation for vorticity, Eq. (9), 
results in a less strong result. Namely, the enstrophy 
(mean squared vorticity) turned out to be an adiabatic 
invariant for the transitional turbulence. 

Adiabatic invariance and scaling transitional regime. 
According to the theoretical predictions if one 

starts from a self-similar solution of the equation (1) 
for weak-nonlinear conditions, then the so-called coher- 
ent regime 0,070; '^ill be developed at a certain stage 
of evolution. The first stage of evolution leads to an 
explosive increase of occupation numbers in the long- 
wavelength region of wavenumber space where the order- 
ing process takes place. From the beginning of the coher- 
ent stage of the evolution we will call the long-wavelength 
region of wavenumber space as pre-condensate fraction, 
while the rest of the wavenumber space we will call as 
above-condensate fraction. 



At a certain time, close to the blow-up time of the self- 
similar solution, the coherent regime sets in. After this 
time the system has a certain transitional turbulent pe- 
riod. At the end of this transitional period appearance of 
a well-defined tangle of quantized vortices indicates the 
final (condensate) stage of the evolution. Therefore, gen- 
eralized vorticity exchange between the pre-condensate 
and above-condensate fractions in the transitional period 
seems to be crucial for the process of the condensate for- 
mation in the KS-scenario. 

In the case of space isotropy it is convenient to deal 
with an angle-averaged occupation numbers spectrum Nk 
in variable k = |k| in the wavenumber space 



N 



Nkdk 



(10) 



The total number of particles associated with the pre- 
condensate fraction is TV^'^ — Jf.,^^. Nl.dk', where kc 
is the wavenumber scale separating (approximately) the 
pre- and above-condensate fractions (for its determina- 
tion from numerical simulations see next section). Then 
the total number of particles associated with the above- 
condensate fraction is: N"-'' = N — N^'^. From the dy- 
namical conservation law for the total number of parti- 
cles, N = const, we obtain 



dNP'' 
dt 



dN" 
dt 



(11) 



Then we can define a single exchange rate of the particles 
number, £„ as 



(12) 



Let us also consider angle-averaged enstrophy spec- 
trum IIa; 









dNP" 




dt 




dt 



nkdk 



(13) 



Then, the enstrophy associated with the pre-condensate 
fraction is fl^'^ = Jf.,^^. n^dfc', whereas enstrophy as- 
sociated with the above-condensate fraction is: fi"'^ = 
(w^) — riP'^. The intensive particle flux to the long- 
wavelength region of wavenumber space (where the pre- 
condensate is formed) involves the enstrophy one. Unlike 
the total number of the particles, the average enstrophy 
{lu'^) is not an exact invariant of the motion. Though, 
characteristic time scale of the enstrophy exchange be- 
tween the pre- and above- condensate fractions is ex- 
pected to be much smaller than the characteristic time 
scale of the (w^) evolution Therefore, the (w^) 

can still be considered (approximately) as an 'adiabatic 
integral' for the exchange process. As it is used for the 
adiabatic processes this statement can be formalized as 
following: 



dvi<'<'/dt\ d\{uj'^)/dty \dVLP^/dt\ \d{uj^)/dt 



(14) 
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(actually, only one of the inequalities (14) is sufficient 
for further consideration). For the nonlinear system the 
principal problem is invariance of the characteristic time- 
scales in (14) to the rescaling ip Xip (otherwise the in- 
equalities (14) are meaningless). It can be readily shown 
that the cancellation of the nonlinear terms in the dy- 
namical equations (8), (9) provides such invariance for the 
enstrophy based time-scales in Eq. (14). 

Since for the transitional turbulence Jl"'^, ^l^'^, and (lj^) 
are still of the same order, then we obtain from Eq. (14) 
following inequalities for the rates 



(15) 





< 






d{LJ^) 


< 


dnp" 


dt 


dt 




dt 


dt 



From n""" = (tj2) _ n'P" and Eq. (15) we obtain 



dnP'' 
dt 



dfl" 



dt 



(16) 



and now we can define a single exchange rate of the en- 
strophy, in full analogy with Eq. (12) 
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FIG. 1: The cumulative number of particles Nk ~ ^k'<k 
against Infc. The data are taken from d for t = 600 (the 
beginning of the transitional stage) . The straight line is drawn 
in order to indicate agreement with Eq. (19). 







dQP" 


dt 




dt 



(17) 



i.e. the adiabatic invariance of the enstrophy substitutes 
its exact invariance in the case of the transitional turbu- 
lence. 

For the weak nonlinearity one can expect that the 
angle-averaged spectrum Nk for sufficiently large k is pro- 
portional to En (cf. [1^,0 )• This can be also valid for 
the transitional turbulence for sufficiently large k. Then 
with the additional dimensional parameter (which has 
dimension T^^) one obtains scaling law for sufficiently 
large k from the dimensional considerations 



Nk 



-1/2 k- 



(18) 



External dissipation or forcing can make the esti- 
mat es ( 14) invalid and, hence, destroy this scenario (cf 
[1H)[12])- Though, in the case of a sufficiently weak and 
linear external dissipation (forcing) the above considera- 
tion can be still valid. In this case the exact invariance 
of the total number of particles N can be replaced by 
its adiabatic invariance (similarly to the enstrophy) and 
the inequalities similar to the Eq.(14) can be made. 
Due to linearity of the dissipation (forcing) these in- 
equalities will be still invariant to the rescaling: ip ^ Xtp. 

Comparison with numerical simulations. In paper [ij 
a large scale </iree-dimensional numerical simulations of 
the Gross-Pitaevskii equation (1) was performed in order 
to reveal all three stages of the evolution from weak tur- 
bulence to superfluid turbulence (the Bose-Einstein con- 
densate) with a tangle of quantized vortices. While the 
observed evolution for t < 600 exhibits the well defined 
self-similar weak turbulence, the period 600 < t < 1000 
was identified in pj with the transitional turbulence. In 
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FIG. 2: As in Fig. 1 but for t = 1000 (the end of the transi- 
tional stage). 



Figs. 1 and 2 one can see the cumulative number of par- 
ticles Nk = Tik'<:k N'l^, which shows how many particles 
have momenta not exceeding k at t — 600 (the beginning 
of the transitional stage, Fig. 1) and at t — 1000 (the 
end of the transitional stage. Fig. 2). 

Using scaling (18) one can estimate this number for 
sufficiently large k as 



Nk^^ 



k'<k 



Nk 



^k^<k'<k ik')-' - Hk/kc) (19) 



where kc is the wavenumber scale separating (approxi- 
mately) the pre- and above-condensate fractions. We use 
the semi-log scales and the straight lines in the Figs. 1,2 
in order to indicate agreement with the Eq. (19). After 
the formation of the quasi-condensate at t > 1000 (third 
stage), the distribution of particles acquires a bimodal 
shape P|. The value kc can be defined at crossing point 
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of the straight hne in Figs. 1,2 (indicating Eq. (19)) and 
the horizontal axis. 

Conclusions. In the above studied KS-scenario the 
particle flux to the long-wavelength region of wavenum- 
ber space (where the pre-condensate is formed) initiates 
and carries the generalized vorticity one. At the transi- 
tion turbulence stage the intensive exchange (transfer) of 
the particles (generalized vorticity) between the pre- and 
above-condensate fractions of the wavenumber space re- 
sults in an approximate adiabatic invariance of the gen- 
eralized enstrophy (the average enstrophy dynamics is 
much slower than the exchange dynamics). On the one 
hand, this adiabatic invariance together with exact con- 



servation law for the total number of particles results in 
a scaling law for the particles number spectrum Eq. (18), 
on the other hand it leads to transition from completely 
disordered situation at the initial weak turbulence to lo- 
cal superfluid order with the tangles of quantized vortices 
(condensate). The scaling law (18) is a quantitative ev- 
idence of this scenario. In this scaling law the initial 
weak nonlinearity is driven by the generalized vorticity 
exchange (under the enstrophy adiabatic invariance) to 
the condensate state that involves significantly different 
length scales and results in the scaling on the transitional 
turbulent stage of the condensate formation. 



[1] N. G. BerlofTl and B. V. Svistunov, Phys. Rev. A 66, 
013603 (2002). 

[2] Yu. Kagan and B.V. Svistunov, Phys. Rev. Lett. 79, 3331 
(1997). 

[3] B.V. Svistunov, J. Moscow Phys. Soc. 1, 373 (1991). 
[4] Yu. Kagan, B.V. Svistunov, and G.V. Shlyapnikov, Zh. 

Eksp. Teor. Fiz. 101, 528 (1992) [Sov. Phys. JETP 75, 

387 (1992)]. 

[5] Yu. Kagan and B.V. Svistunov, Zh. Eksp. Theor. Fiz. 

105, 353 (1994) [Sov. Phys. JETP 78, 187 (1994)]. 
[6] B.V. Svistunov, arXiv:cond-mat/0009368 (2000) 
[7] E. Levich and V. Yakhot,^. PhysrATMath. Gen. 11, 

2237 (1978). 

[8] V.E. Zakharov and S.V. Nazarenko, Physica D, 201 203 
(2005). 

[9] Y. V. Lvov, S. Nazarenko, and R. West, Physica D 184, 
333 (2003). 

[10] N.B. Kopnin, Theory of nonequihbrium superconductiv- 
ity (Clarendon Press, Oxford, 2001). 

[11] M.J. Davis, S.A. Morgan, and K. Burnett, Phys. Rev. 
Lett. 87, 160402 (2001). 

[12] M. Kobayashi and M. Tsubota, Phys. Rev. Lett. 94, 
065302 (2005). 



[13] T. Araki, M. Tsubota, and S.K. Nemirovskii, Phys. Rev. 

Lett. 89, 145301 (2002). 
[14] V. L. Ginzburg and L. P. Pitaevskii, Sov. Phys. JETP 

7, 858 (1958); L. P. Pitaevskii, Sov. Phys. JETP 13, 451 

(1961); E. P. Gross J. Math. Phys. 4, 195 (1963). 
[15] D. V. Semikoz and I. I. Tkachev, Phys. Rev. D, 55, 489 

(1997). 

[16] C.F. Barenghi, R.J. Donnelly, and W.F. Vinen (Eds.), 
Quantized Vortex Dynamics and Superfluid Turbulence 
(Lecture Notes in Physics, Vol. 571), Springer- Verlag, 
2001. 

[17] W.F. Vinen, and J.J. Niemela, Quantum turbulence. J. 
Low Temp. Phys. 128, 167 (2002). 

[18] G. P. Bewley, D. P. Lathrop, and K. R. Sreenivasan, Na- 
ture 441, 588 (2006). 

[19] A. Bershadskii, J. Stat. Phys., 128, 721 (2007). 

[20] J-P. Laval, B. DubruUe and S. Nazarenko, Phys. Fluids, 
13, 1995 (2001). 

[21] A. Dyachenko, and G. Falkovich, Phys. Rev. E 54, 5095 
(1996). 

[22] S. Nazarenko and M. Onorato, Physica D, 219, 1 (2006); 
J. Low Temp Phys. 146 31 (2007). 



